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Ch ' Abstract 

^^ ! Numerical investigation of critical exponents on a hypercubic lattice with L'^ 

random sites with L up to 33 and c? up to 7 showed that above the critical dimension 
the phase transitions in Ising model and percolation are not alike. 
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The upper critical dimension for percolation is 6, instead of 4 for the Ising model. 
Percolation and Ising models are special cases of the Potts model, and the five- 
(^ ■ dimensional Ising case corresponds roughly to percolation in seven dimensions [Ij. 

^ . In the usual Ising model Cheon et al.0 based on a theory of Chen and Dohm |3] 

\^ ', calculated the spontaneous equilibrium magnetisation for d = 5. They used finite- 

O I size scaling in the variables L^ and t = {T — Tc)/Tc according to Ref. 3, the t — l/L"^ — 

S ' plane is similar to that used in fig. 1 of Ref. 3. They found different values of the ex- 

ponent for the spontaneous magnetization in different directions by approaching the 
critical point 1/L = 0, t = 0. Analogous to their previous findings, we study numer- 
ically for random site percolation using the standard Hoshen-Kopelman algorithmic 
'^ \ the behavior of the largest cluster size and the second moment of the cluster size 

distribution for d = 7 as well as d = 3 at, below and above Pc- Above Pc means 
O ■ P — Pc = 1/ L^^" and below pc, P — Pc = —l/L^^"^ where the connectivity length 

exponent z/ is 0.875 (for d = 3) and 1/2 (for d = 7). 
K> '. The infinite cluster normalized to the sample volume is defined as the probability 

H i -Poo for each site to belong that infinite cluster. This probability behaves as Pqo oc 

(p ~ Pc)^- For fine size L, Pqo oc L^""^. The second moment of the cluster size 
distribution x is defined for infinite networks as x oc (p — Pc)~^, and x oc (p — Pc)^ 
for finite size. Using Poo for d = 7 their exponents are summarized for different 
directions as follows: 

(1/L2 = 0, p>p,);(3^1,^ = l (la) 

(1/L2 > 0, p = pc) -,0 = 4.2,0' = 1.8 (16) 

{p-p, = l/L^);D = 5,D' = 2 (Ic) 
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{p-p, = -l/L^)-D = A.2,D' = 2 
and using the relations |1]: 

D = d-(3/v 
D' = -f/iy 



(Id) 

(2) 
(3) 



we obtain the same exponents in the directions a and c. For the Ising model tested 
by Ref. 2, the powers for the directions 6, c and d are respectively 1/L^/'', 1/L and 

1/L3/2. 

Now let us describe the behavior of these exponents in detail. We plot in figure 
1 the maximum cluster exponent D for d = 7 along directions 6, c and d. This 
exponent seems to be close to 4.2 at Pc and also below this critical point with very 
large fluctuations, whereas above this point the exponent seems to be very close to 5 
with very small fluctuations. The fluctuations can be explained by the existence of 
a very large number of small clusters below pc and the critical behavior at this phase 
transition, while above Pc the largest cluster size is very large so that it fluctuates 
only slightly. 







Largest cluster at pc(X), above pc (squares) and below pc (circles) 
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Figure 1: Semi-log plot of the maximum cluster scaling exponent D 
for d — 7. We use sizes from 7 until 33, 64 x 10^ samples for L — 7, 
64000 samples for L = 13, 12800 samples for L ^ 17, 6400 samples for 
L = 19, 60 samples for L = 25 and 10 samples for L — 33. 

We show in figure 2 the exponent D' for the second moment of the cluster size 
distribution (ignoring the largest cluster). At both sides of Pc this exponent is close 
to 2, while it seems to be close to 1.8 at Pc with large fluctuations for the same 
reasons as in figure 1. We cannot compare our data with since there all spanning 







second moment of the cluster size distribution at pc(X), 
above pc (squares) and below pc (circles) 
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Figure 2: Semi-log plot of the exponent D' for d — 7 and the same sizes and 
statistic cited in figure 1 

clusters, not one largest cluster, is studied. 

In order to check the theoretical prediction of D and D' for dimension 3 using 
our Monte Carlo simulation we calculated with less computational efforts these 
exponents for L up to 401 and 500 samples. They seem to not depend on the 
direction and fluctuate strongly around 2.5 and 2 not far from the theoretical values 
2.53 and 2.06. 

Finally, we check the theoretical values of P and 7 along the direction a. We 
show in flgures 3 and 4 the exponent P for the largest cluster and the exponent 7 
for the second moment of the cluster size distribution for Pc approached from above 
for (i = 7 at fixed large L. In figure 3 the exponents cannot be distinghuished above 
p — Pc = 0.008 for all sizes used. But close to pc they differ for /3 by several percent 
from each other and by about 20 percent from the theoretical value 1. In contrast, 
in figure 4, the exponent 7 is in a good agreement with the theoretical value 1 close 
to Pc. 



Conclusion: 

We showed in this work that percolation and Ising models behave differently 
above the critical dimensions in the size range investigated. Because of finite-size 
effects, the exponents should not be expected to agree perfectly with theory. 
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Largest cluster for L=25(+),28(X),29(squares),31 (circles),33(up triangies) 
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Figure 3: The largest cluster exponent (3 of the infinite cluster versus p — Pc 
ioT d = 7. L ~ 33 is the largest size we can simulate, 25 samples for L ~ 25 
and one sample for the other sizes. 



second moment of tfie cluster size distribution 
for L=25(+),28(X),29(squares),31(circies),33(up triangles) 
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Figure 4: The second moment of the cluster size distribution exponent 7 
versus p — Pc for d = 7 with the same sizes and statistic of figure 3 
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